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摘 要 : 青海 省 属于 全 国 四 大 牧区 之 一 ,及 时 监测 草地 植被 长 势 、 准 确 估算 牧草 产量 对 青海 牧区 可 
持续 发 展 与 生态 保护 具有 重要 意义 。 草 地 产 草 量 遥 感 估算 主要 基于 植被 指数 与 地 面 实测 数据 的 统 
计 关 系 ,但 是 估算 涉及 植被 指数 、 统 计 模 型 和 建 模 指标 等 因素 ,不 同 组 合 建立 的 估算 模型 的 精度 不 


同 。 本 文 基于 青海 省 MODIS 数据 与 地 面 实测 产 


量 数 据 ,选择 了 6 种 植被 指数 (NDVI、EVI、RVI、 


DVI RDVI MSAVI) \5 种 统计 模型 ( 简单 线性 模型 、 二 次 多 项 式 模型 、 才 函数 模型 、 指 数 函数 模型 、 对 
数 函 数 模型 ) 以 及 3 种 建 模 指标 ( 植被 指数 年 度 最 大 值 VL 、 植 被 指数 生长 季 累 积 值 VE cason-cum > 1 


被 指数 年 度 累积 值 Vwaowm) ,研究 不 同 组 合 下 估算 模型 的 精度 差异 ,并 从 中 选 出 最 优 产 草 量 估算 


模型 ,用 于 估算 青海 省 2015 年 和 2016 年 的 产 草 量 。 结 果 表 明 :(1)6 种 植被 指数 中 ,基于 NDVI 的 
草 量 估算 精度 最 高 ; 非 线性 模型 的 估算 精度 高 于 线性 模型 ,尤其 是 指数 模型 ,适用 于 大 多 数 草地 
类 型 产 草 量 的 估算 ;基于 NDVI 年 度 最 大 值 的 估算 模型 对 大 多 数 草 地 类 型 都 具有 最 高 的 决定 系数 
(R?)。(2) 从 干 重 来 看 ,高 产 草 量 区 ( >1 200 kg * hm“) 主 要 位 于 青海 东部 的 高 寒 草原 ,中 等 产 草 
区 (600 ~1 200 kg > hm“) 位 于 青海 南部 和 东部 的 高 寒 草 原 和 禾 草 草原 ,低产 草 量 区 ( <600 kg 
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“ hm 一) 位 于 青海 西部 和 北部 的 高 寒 草 旬 \ 高 寒 草原 高 寒 荒 漠 和 盐 生 草 旬 。(3) 与 2015 年 相 比 ， 
2016 年 青海 省 干草 总 产量 减少 31.60 x10”t, 减 幅 为 1.36% 。 其 中 , 禾 草 草原 和 高 寒 草 甸 的 减产 幅 
度 最 大 ,而 荒漠 草原 各 盐 生 草 旬 的 产量 则 有 所 增加 。 本 文 可 为 草地 产 草 量 遇 感 估算 的 研究 和 实践 


提供 参考 。 
关键 词 : 
文章 编号 


青海 省 ; 产 草 量 ; BR; MODIS; 植被 指数 


草地 长 势 监 测 、 产 草 量 估 算 和 和 载 冀 能 力 评估 对 
牧区 的 生态 保护 和 可 持续 发 展 具有 重要 意义 "|。 
产 草 量 不 仅 能 直接 反映 草原 植被 生产 力 ,也 是 开展 
草原 载 畜 能 力 评 估 的 基础 ” ” 。 传 统 的 产 草 量 测定 
方法 耗 时 耗 力 ,难以 在 宏观 斥 度 上 大 规模 应 用 。 自 
20 世纪 90 年 代 以 来 ,大 量 研 究 通过 建立 遥感 植被 
指数 和 地 面 样 方 产 草 量 的 统计 模型 ,对 草地 产 草 量 
及 其 区 域 变化 进行 估算 分 析 “-”i。 研 究 表明 , 产 
草 量 与 光谱 植被 指数 之 间 有 很 强 的 关联 性 。 在 各 种 
植被 指数 中 , 归 一 化 差 值 植 被 指数 (Normalized 
Difference Vegetation Index , NDVI) 最 为 常用 ,是 监测 
植被 生长 活动 的 良好 指标 "-" ,也 是 适用 于 低 生物 
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量 环境 的 植被 指数 。 此 外 ,许多 其 它 植被 指数 ,如 比 
值 植被 指数 (Ratio Vegetation Index, RV/) u 重 归 
一 化 差 值 植被 指数 (Renormalized Difference Vegeta- 
tion Index, RDVI) 增强 型 植被 指数 (Enhanced Vege- 
tation Index, EVI) (12, 17-19) 差 值 植被 指数 ( Difference 
Vegetation Index, DVI) .垂直 植被 指数 (Perpendicular 
Vegetation Index, PVI) !®! 修正 型 土壤 调整 植被 指 
žit ( Modified Soil Adjusted Vegetation Index, MSA- 
VI)" 和 优化 土壤 调整 植被 指数 (Optimized Soil Ad- 
justed Vegetation Index , OSAVI) 上 9 等 ,也 都 被 用 来 估 
算 产 草 量 或 生物 量 。 

但 是 基于 统计 关系 建立 的 产 草 量 估算 模型 具有 
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很 强 的 地 域 性 "1, 由 于 自然 条 件 的 差异 ,不 同 的 植 
被 指数 可 能 适用 于 不 同 的 草地 类 型 。 如 LI 等 ' 通 
过 比较 发 现 ,对 于 低地 草 匈 、 山 地 荡 漠 草原 和 山地 草 
fa), RVI 与 鲜 草 产量 的 相关 性 比 NDVI 更 好 ,但 对 于 
平原 荒漠 草原 则 相反 。 而 YU 等 "和 傅 新 宇 等 '" 
的 研究 表明 ,NDVI 和 EV 在 估算 草地 产 草 量 时 的 精 
度 存 在 差异 ,对 于 大 部 分 高 寒 草 地 而 言 ,基于 NDVI 
的 佑 算 精度 更 高 。 类 似 地 , 杨 秀 春 等 也 发 现 , 相 
比 于 其 他 植被 指数 ,NDVI 和 SAVI 更 适合 用 于 监测 


测 产量 数据 和 MODIS 数据 ,选择 了 6 种 植被 指数 
(NDVI EVI RVI DVI, RDVI,MSAVI) 5 种 统计 模型 
(简单 线性 模型 二 次 多 项 式 模 型 . 宕 函数 模型 . 指 
数 函 数 模 型 对 数 函 数 模型 ) 以 及 3 种 建 模 指标 (年 
度 植被 指数 最 大 值 VI... TER ts BE KE RTR 

VI ,casonoun 、 植 被 指数 年 度 累 积 值 Vw, ) 来 研究 不 
同 估算 模型 的 精度 差异 。 首 先 研 究 哪 种 植被 指数 和 
估算 模型 最 适合 用 于 估算 青海 省 不 同 草地 类 型 的 产 
草 量 , 然 后 比较 基于 3 种 建 模 指 标的 产 草 量 估算 精 


北方 农 牧 交 错 带 的 草原 生长 状况 。 此 外 ,在 产 草 量 
估算 的 统计 模型 方面 ,有 研究 表明 非 线性 模型 往往 
比 线性 模型 的 拟 合 效 果 更 好 ,尤其 是 指数 模型 和 宕 
函数 模型 "””。 并 且 前 人 研究 多 直接 采用 植 
被 指数 最 大 值 来 估算 产 草 量 , 关 于 植被 指数 累积 
或 均值 是 否 能 用 于 产 草 量 佑 算 的 研究 相对 较 
18, 23] 

草地 产 草 量 遥 感 估 算 涉 及 植被 指数 、 统 计 模 型 
和 建 模 指标 等 因素 ,不 同 组 合 建立 的 估算 模型 的 精 
度 不 同 。 尽 管 前 人 已 在 此 方面 进行 了 大 量 的 研究 ， 
但 对 于 青海 草原 而 言 , 哪 种 植被 指数 、 哪 种 统计 模 
型 . 哪 种 指标 最 适宜 用 于 产 草 量 估算 还 有 待 进一步 
研究 。 针 对 以 上 问题 ,本 文 基于 青海 省 草地 地 面 实 
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度 差 异 ,最 终 选 出 最 优 产 草 量 估算 模型 用 于 估算 
2015 年 及 2016 年 的 青海 草原 产 草 量 并 进行 变化 分 
析 ,以 期 为 草地 产 草 量 遥 感 估 算 的 研究 和 实践 提供 
参考 。 


1 研究 数据 与 方法 


1.1 研究 区 概况 

青海 省 位 于 中 国 西北 内 陆 青藏 高 原 东 北部 ,全 
省 下 辖 2 个 市 6 个 州 (图 1), 地 理 坐 标 介 于 
89°35'E ~ 103°04’E 31°39'N ~39°19'N 之 间 。 全 省 
平均 海拔 在 3 000 m 以 上 , 属 高 原 大 陆 性 气候 ,夏季 
短暂 ,年 均 温 在 -6~9 人 之 间 , 平 均 年 降水 量 在 
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1 青海 省 草地 类 型 及 地 面 样 点 分 布 图 


Fig.1 Grassland types and the ground survey samples in Qinghai Province 


Ba ESE: FET MODIS 的 青海 草地 产 草 量变 化 遥感 分 析 


250 ~550 mm 之 间 。 青 海 省 是 我 国 主要 牧区 之 一 ， 
天 然 草地 面积 占 我 国 草地 总 面积 的 10.72% , 达 
4.19 x10’ hm 。 其 中 ,可 利用 草地 面积 约 为 3.87 x 
10 hm? ,主要 分 布 在 东北 部 的 祁连山 地 西南 部 的 
青 南 高 原 以 及 西北 部 的 柴 达 木 盆 地 区 , 主要 草地 类 
型 包括 高 寒 草原 、 禾 草草 原 、 高 寒 草 铅 、 高 寒 荒 漠 
2 

1.2 数据 来 源 

1.2.1 地 面 调查 数据 ”人 研究 区 2015 年 及 2016 年 
的 产 草 量 地 面 实测 数据 来 源 于 青海 省 草原 监理 中 心 
的 地 面 调查 样 方 ,所 有 样 方 的 调查 时 间 集 中 在 每 年 
的 7 月 20 日 ~8 月 31 日 之 间 , 该 时 期 是 草地 植被 
生长 最 旺盛 的 时 期 。 调 查 时 ,矮小 灌木 及 草本 植物 
的 样 方 大 小 为 1 mxl m, 灌 木 及 高 大 草本 样 方 的 大 
小 为 10 m x10 ml 。 调 查 内 容 包括 经 纬度 .草地 
类 型 . 产 草 量 鲜 重 和 样 方 平均 海拔 高 度 等 。 草 地 类 


型 图 (图 1) 源 于 2007 年 出 版 的 《中 华人 民 共 和 国 植 
被 图 (1:1 000 000)》。 


为 保证 地 面 实测 数据 的 可 靠 性 ,本 文 首先 将 记 
录 不 完整 的 样 方 标记 为 无 效 样 方 ,然后 再 根据 样 方 
产 草 量 的 总 体 标准 差 ,将 出 现 较 大 偏差 的 样 方 标记 
为 异常 样 方 , 将 两 者 一 并 剔除 。 符 出 现 2 个 或 多 个 
样 方 落 在 同一 像 元 内 的 情况 , 则 取 其 平均 值 作为 一 
个 样 点 使 用 。 经 过 筛选 最 终 得 到 2015 年 的 有 效 样 
点 数 为 234 个 ,2016 年 的 有 效 样 点 数 为 352 个 , 具 
体 见 表 1 和 图 1。 

因 沼泽 草 旬 没有 地 面 实测 数据 ( 表 1) , 且 其 面 
积 所 占 比 例 非 常 小 (图 1) ,所 以 将 其 排除 在 后 续 分 
析 中 。 尽 管 高 寒 项 漠 缺 少 地 面 实测 数据 ,但 其 所 占 
面积 比较 大 ,由 于 它 与 高 寒 草原 的 地 理 位 置 和 植被 
特征 接近 ,因此 采用 高 寒 草原 的 估 产 模型 估算 其 产 


ya. 
FALE 


表 1 各 草地 类 型 地 面 样 点 数量 
Tab.1 Number of field samples for each grassland type 


草地 类 型 2015 年 样 点 数 2016 年 样 点 数 
高 寒 荒漠 0 0 
禾 草 草原 28 54 
荒漠 草原 7 20 
高 寒 草 原 59 74 
盐 生 草 多 5 8 
高 寒 草 向 135 196 
沼泽 草 多 0 0 

总 计 234 352 


1.2.2 遥感 数据 本 文 结合 了 2 种 MODIS 的 产品 
数据 - MOD09Q1 和 MOD09A1。MOD09Q1 和 
MODO9AL 的 时 间 分 辨 率 均 为 8 d, 但 前 者 只 包含 红 
光波 段 和 近 红 外 波段 的 反射 率 波 段 ,空间 分 辩 率 为 
250 m ,而 后 者 包含 蓝光 、 绿 光波 段 在 内 的 一 共 7 个 
波段 的 反射 率 波段 ,空间 分 辩 率 为 S00 m。 本 文 将 
从 MODOOAL 数据 中 提取 的 蓝光 波段 重 采 样 至 250 
m, 然 后 结合 MOD09Q1 数据 中 的 红 光 波段 和 近 红 外 
波段 ,计算 了 6 种 植被 指数 (NDJ EVI, RVI, DVI, 
RDVI MSAVI) ,计算 公式 如 表 2 所 示 。 


表 2 6 种 植被 指数 计算 公式 


Tab.2 Equations of the six vegetation indices 


植被 指数 公式 
NDTTL26] NDVI -PME 一 AR 
PNIR 十 DR 
] 2. SPwR Pr) 
Eyl? EVI = 
Prin +6. Opr —7. Spr +1 
Rv) RVI = 各 性 
PR 
bys! DVI = pyr -Prl 28] 
RDT7L29] RDVI = PAR PR 
Pur 十 DR 


MSAVI'*0] MSAVI = 本 | 2pwR +1- VM (2pnie +1) -8(pwn -Pr) | 


注 :pwr PR 和 pr 分 别 表示 近 红 外 波段 、 红 光波 段 和 蓝光 波段 的 反 
HK 


受 海拔 高 度 和 水 汽 的 影响 ,研究 区 的 大 气 状况 
多 变 ,数据 中 存在 的 噪声 会 影响 产 草 量 估算 的 准确 
性 。 因 此 ,本 文通 过 线性 插值 对 突然 上 升 或 下 降 的 
点 进行 替换 ,然后 计算 了 以 下 三 项 建 模 指 标 : (1) 全 
年 植被 指数 最 大 值 , 记 为 W,,.。(2) 生 长 季 植 被 指 
数 累 积 值 , 记 为 VLeasoncm i (3 ) 全 年 植被 指数 累积 
值 , 记 为 Vin。 青 海 省 草地 一 般 在 每 年 的 5 月 
上 名 返青 ,在 每 年 的 9 月 底 进 入 枯黄 期 ,因此 本 文 将 
草地 的 生长 季 定 义 为 5 ~9 月 。 

1.3 研究 方法 

1.3.1 产 草 量 估算 方法 首先 按 不 同 草 地 类 型 , 建 
立 6 种 植被 指数 (NDVI、EVI、RVI、DVI、RDVI、MS4- 
VI) 的 WV, 与 对 应 地 面 样 点 鲜 草 产量 的 回归 模型 ， 
从 中 选 出 决定 系数 (R ) 最 高 的 植被 指数 作为 最 优 
植被 指数 。 然 后 ,将 最 优 植被 指数 的 三 项 建 模 指 标 
(VL + VE season-cum Tem) 和 5 种 统计 回归 模型 ( 线 
性 模型 .多项式 模型 AC PR I FF PRI A AT 
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RI 各 草地 类 型 干 鲜 比 折算 系数 
Tab.3 Conversion coefficient for fresh grass yield 


to hay yield for each grassland type 


草地 类 型 折算 系数 
CESARA 1.0/2.5 
禾 草草 原 1.0/3.0 
荒漠 草原 1.0/3.0 
高 寒 草 原 1.0/3.0 
Eh AE BY fy 1.0/3.2 
AXHA 1.0/3.2 


数 函 数 模 型 ) 两 两 组 合 ,建立 鲜 草 产量 估算 模型 ,再 
采用 回归 方程 决定 系数 (R ) RMSE 和 REE 来 评价 
各 模型 精度 ,从 中 选 出 最 优 模型 估算 鲜 草 产量 。 最 
后 ,根据 《中 国 草地 资源 》” 中 的 干 鲜 比 折算 系数 
( 表 3) ,将 鲜 草 产量 换算 成 干草 产量 并 进行 分 析 。 
1.3.2 模型 验证 在 建立 产 草 量 估算 模型 前 ,将 样 
点 数据 分 为 建 模样 本 和 检验 样本 两 部 分 ,其 中 建 模 
样本 约 占 4/5 ,检验 样本 约 占 1/5, 各 草地 类 型 建 模 
样 点 数 及 检验 样 点 数 具体 如 表 4 所 示 。 其 中 ,有 2 
种 草地 类 型 由 于 样 点 数量 过 少 , 无 法 对 其 模型 精度 
进行 检验 ( 表 4) 。 在 验证 模型 精度 时 ,采用 回归 方 
FEW thee AMR’ .RMSEI2 23 和 REEMA 5 三 项 
指标 来 评价 模型 精度 ,公式 如 (1) 和 (2) 所 示 : 


m 


式 中 :了 为 第 i 个 样本 的 实测 鲜 草 产量 (kg + hm?) 5 
Y, 为 第 i 个 样本 的 估算 鲜 草 产量 (kg + hm-?); 了 为 
Y, 的 平均 数 ;n 为 检验 样本 数 。 

表 4 各 草地 类 型 建 模 及 检验 样 点 数 


Tab.4 Modeling and validation samples of each 
grassland type 


建 模样 点 数 检验 样 点 数 总 样 点 数 
草地 类 型 2015 年 2016 年 2015 年 2016 年 2015 年 2016 年 
禾 草草 原 ”22 43 6 11 28 54 
荒漠 草原 7 16 O# 4 7 20 
高 寒 草 原 48 58 11 16 59 74 
盐 生 草 旬 5 8 O# O# 5 8 
Beefy 109 155 26 41 135 196 
总 计 191 280 43 72 234 352 
不 预 留 


注 :# 表 示 由 于 没有 足够 的 样 点 , 故 将 所 有 样 点 用 于 建 模 使 
检验 样 点 


2 结果 与 分 析 


2.1 最 优 产 草 量 估算 植被 指数 

基于 6 种 植被 指数 全 年 最 大 值 (VL, ) , 本文 建 
立 了 2015 年 和 2016 年 各 草地 类 型 的 鲜 草 产量 估算 
模型 。 以 2016 年 为 例 ,各 估算 模型 的 决定 系数 
(及 ) 如 图 2 所 示 。 在 大 多 数 情况 下 ,基于 NDV gan 
RDVI 和 Rs 的 产 草 量 估算 精度 更 高 。 从 不 同 草 


an , 、 MARR FARGER, JOEA 
= a (1) SRA, NDVI a 都 是 产 草 量 估算 的 最 优 植被 指数 。 
REE = RMSE/Y (2) ”其 中 , 禾 草草 原 的 二 次 多 项 式 模型 R 最 高 ,为 0.79; 
一 6 一 NDVImax —s— EVImax —+— RVImax —s— DVImax —*— RDVImax —*— MSAVImax 
1:00 OREHE ad OFRE 100 (© 高 寒 章 原 
0.85 0.85 0.85 
ay 0-70 = l A a 0-70 0.70 
0.55 | Ee——#—y 0.55 0.55 
0.40 ——— 0.40 0.40 
0.25 0.25 0.25 
M1 M2 M3 M4 M5 M1 M2 M3 M4 M5 M1 M2 M3 M4 M5 
模型 模型 模型 
1.00 OREHA 1.00 (e) 高 寒 草 多 
0.85 a 0.85 
0.70 a 0-70 ES 
0.55 0.55 
0.40 0.40 
0.25 0.25 
M1 M2 M3 M4 M5 M1 M2 M3 M4 M5 
模型 模型 


注 :MI 线性 模型 ;M2 二 次 多 项 式 模型 ;M3 短 函 数 模型 ;M4 指数 函数 模型 ;M5 对 数 函 数 模型 
图 2 2016 年 植被 指数 全 年 最 大 值 与 鲜 草 产量 的 回归 模型 决定 系数 


Fig.2 R? of regression model between the annual maximum of vegetation index and fresh grass yield in 2016 
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—e— NDVI. —e— NDV ow oum > NDVI wa om 
ai (a) 禾 草 草原 1400 (a) 禾 草 草原 (a) RAB 
0.80 E 1100 x 
= oe u N 
=% 0.60 T $p 800 B 70 
省 一 = æ 65 
0.40 2 500 & 60 
0.20 = 200 55 
Ml M2 M3 M M Ml M2 M3 M M Ml M2 M3 M M5 
模 型 模 型 模 型 
opas 1400 b) 荒漠 草原 = (b,) HRH 
0.80 E 1100 zx 
S% 
œ 0.60 We Èp 800 B 70 
0.40 2 500 I 60 
0.20 = 200 55 
Ml M2 M3 M M5 M M2 M3 M M5 Ml M2 M3 M4 M5 
模 型 模 型 模 型 
casam 1400 C) ARER = (C) ARER 
0.80 E 1100 x 
~ 75 
x 0.60 Èp 800 A = 70 RN 
S R 65 
0.40 2 500 sð 
0.20 = 200 55 
Ml M2 M3 M M5 Ml M2 M3 M M5 Ml M2 M3 M M5 
模 型 模 型 模 型 
1.00 1400 85 
(a) BAA ® Cd) BAA i (dy) ARAA 
0.80 E 1100 f 
ce 0.60 — ia ma B 800 | gy x 2 2 70 
3 zi 65 
0.40 2 500 F 60 
0.20 200 55 
Ml M2 M M M Ml M M3 M M Ml M2 M3 M4 M5 
模 型 模 型 模 型 
(C) 盐 生 草 多 
0.80 SS 
ce 0.60 
9.20 图 3 2016 年 基于 不 同 NDV1 指 标的 估算 模型 精度 
0.20 Fig.3 Validations of estimation models based on different metrics in 2016 
M M2 M3 M M 
模 型 


荒漠 草原 的 震 函 数 模型 尺 最 高 ,为 0.81。 对 于 高 
寒 草 原 和 高 寒 草 旬 ,MD 了 .在 大 多 数 情况 下 的 估算 
精度 更 高 。 其 中 ,基于 指数 模型 的 R 最 高 ,分 别 为 
0.85 和 0.72。 对 于 盐 生 草 馈 而 言 ,在 大 多 数 情况 下 
BR DVL,,, 外 ,其 他 植被 指数 的 R 非常 接近 。 其 中 ， 
基于 DV 的 二 次 多 项 式 模型 R 最 高 , 达 0.95。 综 
合 来 看 ,在 大 多 数 情况 下 非 线性 模型 比 线性 模型 的 
拟 合 效果 更 好 ,基于 NDVI 的 估算 模型 相 较 于 其 他 
植被 指数 具有 更 高 的 决定 系数 (及 )o 2015 年 所 得 
结果 的 趋势 与 2016 年 一 致 。 
2.2 最 优 产 草 量 估算 模型 

根据 上 述 分 析 , 进一步 比较 了 基于 NDVI 的 全 


年 最 大 值 ( NDVL,, ) 、 生 长 季 累 积 值 (NDVI ) 和 
4E BE RARI (NDVI mnualcum ) 的 产 草 量 估算 模型 精度 。 
以 2016 年 为 例 ( 图 3) ,对 于 荒漠 草原 高 寒 草 原 .高 
Fe Ee fa) A ERE EE fa), 基于 内。 x Wl ceason-cum > 
Ja 这 3 PAHS Hay BY PA Be Eat rk AARAA R? GEN 
接近 ,而 对 于 禾 草 草原 ,基于 NDVL, 的 产 草 量 估算 
模型 的 R 明显 高 于 基于 另外 2 种 指标 的 估算 模型 。 
此 外 ,在 大 多 数 情况 下 ,基于 NDVL,, 的 估算 模型 的 
RMSE 低 于 基于 NDVI ocason-cum 和 基于 NDVLiwoaom 的 
估算 模型 ,而 模型 精度 (1 - REE) 高 于 基于 NDV ea 
wum 和 基于 NDJam 的 估算 模型 。 因 此 ,综合 
R? RMSE 和 模型 精度 这 三 项 指标 ,基于 NDVI, HE 


tó wy 


佑 算 模 型 效果 最 好 ,最 适合 青海 省 的 牧草 产量 估算 。 
2015 年 所 得 结果 的 趋势 与 2016 年 一 致 。 

通过 上 述 分 析 选 出 2015 年 及 2016 年 各 草地 类 
型 的 最 优 佑 产 模 型 。 考 虑 到 使 用 不 同 植被 指数 和 统 
计 模 型 可 能 会 产生 潜在 的 不 确定 性 ,因此 选择 最 优 
佑 产 模型 时 ,在 精度 相差 不 大 的 情况 下 采用 相同 的 
植被 指数 和 统计 模型 ,最 终结 果 如 表 5 所 示 。 
2.3 产 草 量 空间 分 布 

基于 上 述 各 草地 类 型 的 最 优 估 产 模型 计算 鲜 草 
产量 后 ,再 乘 以 表 1 中 列 出 的 转换 系数 得 到 干草 产 
量 。2015 年 及 2016 年 青海 省 鲜 草 产量 和 干草 产量 
的 空间 分 布 如 图 4 所 示 。 结 合 图 1 可 知 ,不 同 草地 
类 型 及 不 同 地 区 的 产 草 量 存 在 明显 的 差异 。 产 草 量 
的 空间 分 布 格局 主要 与 草地 类 型 海拔、 地形 等 有 
关 。 从 干 重 来 看 , 产 草 量 较 高 的 区 域 ( >1 200 kg- 
hm“) 主 要 分 布 在 东部 高 寒 草原 ;中 等 产量 区 (600 
~1 200 kg + hm ) 分 布 在 南部 和 东部 的 高 寒 草原 
以 及 禾 草 草原 ;低产 量 区 ( <600 kg hm?) 主要 分 
布 在 西部 和 北部 的 高 寒 草 甸 ,高 寒 草 原 ,高寒 荒 漠 以 
及 盐 生 草 多 。 
2.4 产 草 量 变化 分 析 
2.4.1 不 同行 政 区 的 产 草 量变 化 ”2015 年 和 2016 

青海 省 各 州 市 产 草 量 及 其 变化 情况 如 表 6 所 示 。 

以 2016 年 为 例 ,最 大 的 3 个 干草 产 区 是 南部 和 西部 
的 玉树 州 . 果 洛 州 和 海 西 州 ,干草 产量 共计 1.55 x 
107t, 占 全 省 总 产量 的 67.63%。 此 外 ,有 6 个 地 区 
的 平均 干草 单产 超过 500 kg- hm ,主要 位 于 青海 


省 东南 部 ,包括 黄南州 .西宁 市 、 海 北 州 、 海 东 市 、 果 
洛 州 和 海南 州 ;平均 干草 单产 低 于 500 kg + hm -的 
地 区 主要 位 于 西北 部 的 海 西 州 和 玉树 州 。 尽 管 青 
海 省 西部 地 区 的 干草 产量 占 全 省 很 大 比例 ,但 由 
于 地 形 气候 等 自然 条 件 恶 劣 ,其 干草 单产 远 低 于 
东南 部 地 区 ,而 东南 部 地 区 虽然 干草 单产 较 高 ,但 
由 于 草原 面积 相对 较 小 ,其 干草 产量 仅 占 全 省 的 
一 小 部 分 。 

图 5 表示 了 青海 省 2015 一 2016 年 干草 产量 的 
变化 情况 。 与 2015 年 的 产量 相 比 ,2016 年 海 东 市 
和 黄南州 的 产量 显著 增加 ,超过 5%。 如 图 Sb 所 
示 , 这 两 个 地 区 干草 产量 增加 大 于 150 kg . hm 一 的 
像 元 比例 超过 35% 。 此 外 ,海南 州 东 部 、 玉 树 州 南 
部 、 果 洛 州 东部 以 及 海北 州 的 一 些 地 区 也 出 现 了 显 
著 的 增长 。 从 绝对 产量 来 看 ,黄南州 增产 量 最 高 , 达 
13.99 x 10° 1, 而 海 西 州 减产 量 最 高 , 达 43. 98 x 10° 
t。 减 产地 区 主要 集中 在 海 西 州 东部 、 果 洛 州 北部 和 
海南 州 西部 ,其 中 海 西 州 有 75% 以 上 的 草地 产量 下 
降 ,干草 总 产量 下 降 了 12. 98% 。 尽 管 青海 省 部 分 
地 区 的 产量 大 幅 增加 ,但 由 于 海 西 州 及 玉树 州 北部 
的 产量 普遍 下 降 ,导致 2016 年 全 省 干草 总 产量 比 
2015 年 减少 31. 60 x 10* +t, 减产 幅度 为 1.36% 。 
2.4.2 不 同 草地 类 型 的 产 草 量变 化 

青海 省 2015 年 和 2016 年 各 草地 类 型 的 产 草 量 
及 其 变化 情况 如 表 7 和 图 6 所 示 。 以 2016 年 为 例 ， 
前 三 大 产 草草 地 类 型 为 高 寒 草 旬 Te FE HE AA AR BE 
草原 , 共 占 全 省 总 产 草 量 的 95.58% 。 从 2015 年 到 


RS 2015 年 及 2016 年 各 草地 类 型 的 最 优 估 产 模型 
Tab.5 Optimal grass yield estimation models for each grassland types in 2015 and 2016 


年 份 。 ”草地 类 型 。 ”植被 指数 。 ”最 优 模型 模型 表达 式 R ee ee 

2015 年 REEE NDVI pa 二 次 多 项 式  y=2 333.9 -5 800. 9x +10 332x? 0.81 = 541.22 79.05 
荒漠 草原 NDVI 短 函 数 模 型 ” y=3 261. 3x) 8! 0.95 * - - 
高 寒 草 原 NDVI, 指数 模型 y =250.74e3 7 0.83 “* 331.27 80. 24 
盐 生 草 负 NDVI max 指数 模型 y =449.14e3 ** 0.95 ** - = 
AEE fi NDVI. 指数 模型 y =391. 340 0.74 ** 660. 96 71.17 

2016 年 REEE NDVI nax 二 次 多 项 式  y=1 684.7 -4 841x +10 624x? 0.79 = 639.61 76.42 
荒漠 草原 NDVI max TERRORA  y=4 570. 7x! 0.81" 337.93 75.65 
高 寒 草原 NDVI ma 指数 模型 y = 184. 09e4 8 0.85 “* 678. 43 74.72 
盐 生 草 负 NDVI max 指数 模型 少 =594. 22e3 * 0. 81 * - = 
高 寒 草 馈 NDVI max 指数 模型 y = 191. 786° 8* 0.72 “* 658. 08 78.00 

注 :y 表示 鲜 草 产量 / kg + hm? ,x 表示 NDI7, 数 值 范围 为 -1.0 ~ 1.0; xx 址 示 通 过 0.001 显著 性 检验 (P <0.001) ; -表示 无 数据 
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90° E 95° E 100° E 90° E 95° E 100° E 
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90° E 95° E 100° E 90° E 95° E 100° E 
90° E 95° E 100° E 
a Z 
o (d) 20164 o 
S = 
Z z 
3 3 
a Z 
io io 
cm oO 
鲜 草 产量 /kg*hm” saps FHE /kgehm? | 
Z E <500 E 2 000~4 000 Z = <300 Œ 900-1 200 
& | Œ 500~1 000 Æ 4000-6000 0 100 200 400 & | 21300~600 EE 1 200~1 500 0 100200 400 
E 1 000~2 000 E >6 000 So km = 600~900 E >1 500 oS km 
90° E 95° E 100° E 90° E 95° E 100° E 
4 2015 年 及 2016 年 青海 省 鲜 草 产量 和 干草 产量 的 空间 分 布 
Fig.4 Spatial distribution of fresh grass yield and the hay yield in 2015 and 2016 in Qinghai Province 
表 6 2015 年 及 2016 年 青海 省 各 州 市 产 草 量 
Tab.6 Grass yield statistics at prefecture — level in Qinghai Province in 2015 and 2016 
地 区 草地 面积 鲜 重 /10" t 干 重 / 10% t 平均 干草 单产 / kg "hm ” 干 重 绝对 变 变化 幅度 
x 
/10 hm? 2015 4 2016 年 2015 年 2016 年 2015 年 20164 “化 量 /10 t 7 多 
西宁 市 44.15 148.59 152.92 47.94 49.34 1 085.77 1 117.60 1.41 2.94 
海 东 市 76. 14 172.16 192.75 56. 64 63.47 744.34 833. 99 6. 82 12.04 
海北 州 253. 34 714.07 724.37 225.29 228.65 891.76 905. 06 3.36 1.49 
海 西 州 1 417.92 1 051.05 915.99 338.91 294.93 314.71 273.87 -43.98 -12.98 
海南 州 292. 19 665.62 678. 06 214.47 218.83 734. 10 749.02 4.36 2.03 
WH ”1 570.60 2 440.43 2 425.91 768.83 763.69 490. 41 487.13 -5.14 -0.67 
果 洛 州 619. 10 1 603. 18 1 563.92 501.67 489.26 810.79 790. 73 - 12.42 一 2.47 
黄南州 143. 13 531.11 575.67 166. 63 180. 62 1 164.17 1 261.90 13.99 8.39 
青海 省 ”4 416.58 7 326.22 7 229.60 2 320.38 2 288.78 569.90 562.14 -31.60 -1.36 
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图 5 2015 一 2016 年 青海 省 各 州 市 干草 产量 变化 
Fig.5 Hay yield changes at the prefecture-level in Qinghai Province from 2015 to 2016 


R7 2015 年 及 2016 年 青海 省 各 草地 类 型 产 草 量 
Tab.7 Grass yields of different grassland types in Qinghai Province in 2015 and 2016 


Wh ] 类 刑 He Als 

a 2015 年 ”2016 年 ”2015 年 ”2016 年 ”2015 年 ”2016 年 ”2015 年 ”2016 年 《10t % 
高 寒 荒漠 53.69 30.36 26. 15 12.14 10.46 226.88 195.39 0.50 0.44 -1.68 -13.84 
禾 划 草原 219.13 485.55 465.79 161.85 155.26 739.38 709.28 6.71 6.56 -6.59 -4.07 
荒漠 草原 。 55.54 49.90 74.28 16.63 24.76 299.90 446.47 0.69 1.05 8.13 48.87 
高 寒 章 原 1135.15 958.51 892.56 319.50 297.52 282.24 262.82 13.25 12.57 -21.98 -6.88 
盐 生 草 旬 204.11 177.81 222.19 55.57 69.43 278.31 347.76 2.30 2.93 13.87 24.96 
PIM fy] 2749.48 5905.01 5792.19 1845.32 1810.06 671.72 658.89 76.54 76.45 -35.26 -1.91 
总 计 4417.10 7607.14 7473.16 2411.01 2367.50 545.84 535.98 100.00 100.00 -43.51 -1.80 
干草 产量 变化 /kg.hm 量 减 少 。 其 中 , 荒漠 草原 的 干草 产量 增幅 最 大 ,为 


a<-150 e-80~-15000~-80 90~80 a80~150 a>150 


高 寒 荒漠 
ABER 
荒漠 草原 
高 寒 草原 
盐 生 草原 
AREA 


百分比 时 /% 
图 6 2015—2016 年 青海 省 各 草地 类 型 十 草 产量 变化 
Fig.6 Hay yield changes of different grassland types in Qinghai 
Province from 2015 to 2016 


2016 年 ,荒漠 草原 和 盐 生 草 多 的 干草 产量 增加 ,但 
高 寒 荒漠 、 禾 草草 原 .高寒 草 原 和 高 寒 草 向 的 干草 产 


48.87% ,高 涯 荒漠 的 干草 产量 减 幅 最 大 ,为 13. 84% 。 

此 外 ,如 图 6 所 示 ,低产 草地 类 型 (如 鞠 漠 草原 
和 盐 生 草 甸 ) 的 增产 像 元 比例 超过 80% 。 相 比 之 
下 , 禾 草 草原 和 高 寒 草 甸 等 高 产 草地 类 型 的 产 草 量 
略微 下 降 。 其 中 , 禾 草 草原 干草 产量 下 降 超 过 150 
kg + hm ”的 像 元 比例 大 于 25% 。 总 体 来 看 ,青海 省 
干草 产量 表现 出 略微 的 减产 趋势 ,产量 减少 在 0 ~ 
80 kg + hm “的 像 元 占 大 多 数 。 


3 ”结论 与 讨论 


草地 产 草 量 变化 分 析 是 载 冀 能 力 评估 和 草地 气 
候 变化 响应 等 研究 中 非常 关键 的 一 步 。 遥 感 产 草 量 
估算 涉及 到 植被 指数 、 统 计 模 型 、 建 模 指标 、 草 地 类 
型 以 及 不 同 的 处 理 和 建 模 方式 等 要 素 ,这些 要 素 对 
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产 草 量 估算 结果 存在 一 定 的 影响 。 为 了 提高 产 草 量 
估算 精度 ,本 文 基于 地 面 实测 数据 和 MODIS 数据 ， 
选择 了 6 PP AL RIE BX (NDVI, EVI, RVI, DVI, RDVI, 
MSAVI) 5 种 常用 的 统计 模型 (线性 模型 二 次 多 项 
式 模 型 窜 函 数 模 型 指数 函数 模型 对 数 函 数 模 
型 ) 和 3 种 建 模 指标 (Vl, Vl ceason-cum » Vannual-cum ) 进行 
产 草 量 估算 精度 比较 ,在 此 基础 上 选择 最 优 产 草 量 
估算 模型 估算 了 2015 年 和 2016 年 青海 省 草地 的 产 
草 量 ,并 对 产 草 量 的 变化 情况 进行 分 析 。 通 过 比较 
发 现 ,6 种 植被 指数 中 ,基于 NDVI 的 产 草 量 估算 精 
度 最 高 。 这 与 前 人 在 高 寒 草地 进行 研究 所 得 的 结论 
一 致 .…" ,与 其 他 草地 研究 中 得 出 的 结论 不 一 
BNP) ,本 文 进一步 证 实 了 NDVI 确实 是 更 加 适合 
用 来 估算 高 寒 草地 的 产 草 量 的 植被 指数 。 而 非 线性 
模型 的 估算 精度 高 于 线性 模型 ,尤其 是 指数 模型 , 适 
用 于 大 多 数 草 地 类 型 产 草 量 的 估算 。 这 与 前 人 在 不 
同 区 域 不 同 草 地 类 型 进行 研究 所 得 的 结论 非常 一 
Blt 21) ES NDVI 生长 季 累 积 值 和 年 度 累 
积 值 , NDVI 年度 最 大 值 的 估算 模型 对 大 多 数 草地 类 
型 都 具有 最 高 的 决定 系数 (R*) ,说 明 植 被 指数 最 大 
值 对 于 指示 草地 产 草 量具 有 较 好 的 作用 ,能 够 代表 
草地 最 大 的 生产 力 。 

青海 省 草地 产 草 量 的 空间 分 布 格局 主要 与 草地 
类 型 海拔 .地 形 等 因素 有 关 。 干 草 产 量 较 高 的 区 域 
( >1 200 kg + hm?) 主要 分 布 在 东部 高 寒 草 原 ; 中 


等 产 草 量 区 (600 ~ 1 200 kg + hm ) 分 布 在 南部 和 
东部 的 高 寒 草 原 和 不 草草 原 ; 低 产 草 量 区 ( < 600 
kg + hm 一) 主要 分 布 在 西部 和 北部 的 高 寒 草 多 \ 高 
寒 草原 、 高寒 蕊 并 和 盐 生 草 甸 。 与 2015 年 相 比 ， 
2016 年 青海 省 干草 产量 略 有 下 降 , 下 降 了 31.6 x 
10° t, 减 幅 为 1.36% ,产量 增加 的 草地 类 型 主要 为 
充满 草原 和 盐 生 草 和 多 ,而 产量 减少 的 主要 为 条 草草 
原 和 高 寒 草 旬 。 本 研究 不 仅 可 为 青海 牧区 生态 保护 
以 及 可 持续 发 展 提供 参考 ,同时 也 为 产 草 量 冰 感 估 
算 研究 提供 了 有 益 的 补充 。 


参考 文献 ( References ) 


[1] ALII,CAWKWELL F,DWYER E,et al. Satellite remote sensing 
of grasslands: From observation to management [ J]. Journal of 
Plant Ecology ,2016 ,9(6) :649 - 671. 

[2] 韩 其 飞 , 陆 研 , 李 超凡 . 气候 变化 对 中 亚 草地 生态 系统 碳 循环 
的 影响 研究 [ 中. 干旱 区 地 理 ,2018,41 (6):1351 - 1357. [ 
HAN Qifei, LU Yan, LI Chaofan. Impact of climate change on 


[4] 


[5] 


[7] 


[8] 


[11] 


[12] 


[13] 


grassland carbon cycling in Central Asia[ J]. Arid Land Geogra- 
phy ,2018 ,41(6) :1351 - 1357. ] 

王 顺利 , 王 荣 新 ,敬文 茂 ,等 . 祁连山 干旱 山地 草地 生物 量 对 水 
分 条 件 的 响应 [可 .干旱 区 地 理 ,2017,40 (4) :772 - 779. 
[ WANG Shunli, WANG Rongxin, JING Wenmao, et al. Biomass of 


grassland and response to soil moisture on arid mountain land in 
the Qilian Mountains[ J]. Arid Land Geography ,2017 ,40(4) :772 
-779. ] 

XU B, YANG X C,TAO W G,et al. MODIS-based remote sensing 
monitoring of grass production in China[ J]. International Journal 
of Remote Sensing ,2008 ,29(17 -18 ) :5313 - 5327. 

R, DEE. AR AC ER DK et ATT BE BF EAL I]. 
地 理 研 究 ,2009 ,28 (2) :402 - 408. [ XU Bin, YANG Xiuchun. 
Calculation of grass production and balance of livestock carrying 
capacity in rangeland region of northeast China[ J]. Geographical 
Research ,2009 ,28(2) :402 - 408. ] 

徐 剑 波 , 宋 立 生 , 赵 之 重 ,等 . 近 15 a 来 黄河 源 地 区 玛 多 县 草地 
植被 退化 的 遥感 动态 监测 [中 .干旱 区 地 理 ,2012,35 (4) :615 
—622.[ XU Jianbo,SONG Lisheng,ZHAO Zhizhong, et al. Moni- 


toring grassland degradation dynamically at Maduo County in 
source region of Yellow River in past 15 years based on remote 
sensing[ J]. Arid Land Geography ,2012 ,35(4) :615 —622. ] 

LI X W,LI M D,DONG S K,et al. Temporal-spatial changes in e- 
cosystem services and implications for the conservation of alpine 
rangelands on the Qinghai-Tibetan Plateau [ J}. Rangeland Jour- 
nal ,2015 ,37(SI) :31 -43. 

MA W H,FANG J Y, YANG Y H ,et al. Biomass carbon stocks and 
their changes in northern China’s grasslands during 1982 - 2006 
[J]. Science China-Life Sciences ,2010 ,53(7) :841 - 850. 

PIAO S L,FANG J Y,ZHOU L M,et al. Changes in biomass car- 
bon stocks in China’s grasslands between 1982 and 1999 J}. 
Global Biogeochemical Cycles ,2007 ,21 ( GB2002 ) : doi; 10. 1029/ 
2005 GB002634. 

RUSCH G M, ZAPATA P C,CASANOVES F, et al. Determinants 
of grassland primary production in seasonally-dry silvopastoral sys- 
tems in Central America J]. Agroforestry Systems ,2014 ,88 (3) : 
517 - 526. 

WEHLAGE D C,GAMON J A,THAYER D,et al. Interannual var- 
iability in dry mixed-grass prairie yield; A comparison of MODIS, 
SPOT and field measurements [ J]. Remote Sensing, 2016, 8 
(87210) : doi; 10. 3390/rs8 100872. 

YU L,ZHOU L,LIU W,et al. Using remote sensing and GIS tech- 
nologies to estimate grass yield and livestock carrying capacity of 
alpine grasslands in Golog Prefecture, China [ J}. Pedosphere, 
2010 ,20(3) :342 - 351. 

MOHAMMAT A, WANG X H, XU X T,et al. Drought and spring 
cooling induced recent decrease in vegetation growth in Inner Asia 


[J]. Agricultural and Forest Meteorology ,2013 ,178 - 179:21 — 


202006.00132v1 


chinaXiv 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


tó wy 


30. 

PIAO S L,MOHAMMAT A,FANG J Y,et al. NDVI-based increase 
in growth of temperate grasslands and its responses to climate chan- 
ges in China[ J]. Global Environmental Change-human and Policy 
Dimensions ,2006 ,16(4) :340 — 348. 

ZHOU LM, TUCKER C J, KAUFMANN R K,et al. Variations in 
northern vegetation activity inferred from satellite data of vegetation 
index during 1981 to 1999[ J]. Journal of Geophysical Research- 
Atmospheres ,2001 ,106( D17 ) :20069 - 20083. 

LI J L, LIANG T G,CHEN Q G. Estimating grassland yields using 
remote sensing and GIS technologies in China[ J]. New Zealand 
Journal of Agricultural Research , 1998 ,41 (1) :31 -38. 

传 新 宇 , 唐 川 江 , 张 新 路 ,等 . 四 川 草 原 MODIS 数据 的 产 草 量 
估算 [中 .地 球 信息 科学 学 报 ,2013,15(4):611 - 617. [ FU 
Xinyu, TANG Chuanjiang, ZHANG Xinyue, et al. Estimation of 


grass yield based on MODIS data in Sichuan Province, China[ J}. 
Journal of Earth Information Science ,2013 ,15(4) :611 - 617. ] 

杨 淑 起 , 张 文 娟 , 冯 琦 胜 ,等 . 基于 MODIS 逐日 地 表 反 射 率 数 
据 的 青 南 地 区 草地 生长 状况 遥感 监测 研究 [站]. 草 业 学 报 ， 
2016 ,25(8) :14 —26. [ YANG Shuxia, ZHANG Wenjuan, FENG 


Qisheng , et al. Monitoring of grassland herbage accumulation by re- 
mote sensing using MODIS daily surface reflectance data in the 
Qingnan Region[ J]. Acta Prataculturae Sinica ,2016 ,25(8) :14 - 
26. | 

杨 秀 春 , 徐 斌 , 朱 晓 华 , 等 . 北方 农 牧 交 错 带 草原 产 草 量 遥 感 监 
测 模 型 [站 . 地 理 研 究 ,2007,(2) :213 -221.[ YANG Xiuchun, 


XU Bin, ZHU Xiaohua,et al. Models of grass production based on 
remote sensing monitoring in northern agro-grazing ecotone [ J ]. 
Geographical Research ,2007 , (2) :213 -221. ] 

DI BELLA C,FAIVRE R,RUGET F, et al. Remote sensing capa- 
bilities to estimate pasture production in France[ J |. International 
Journal of Remote Sensing ,2004 ,25 (23 ) :5359 — 5372. 
张扬 建 , 范 春 捆 , 黄 珂 ,等 . 遥感 在 生态 系统 生态 学 上 应 用 的 机 
i SPEAR J]. 生态 学 杂志 ,2017,36(3) :809 -823.[ ZHANG 


Yangjian, FAN Chunkun, HUANG Ke, et al. Opportunities and 


challenges in remote sensing applications to ecosystem ecology 
[J]. Chinese Journal of Ecology ,2017 ,36(3) :809 — 823. ] 

罗 玲 , 王 宗明 , 任 春 颖 ,等 . 基于 MODIS 数据 的 松 嫩 草 原 产 草 
量 遥 感 估 算 模型 与 空间 反 演 [J]. 农业 工程 学 报 ,2010 ,26(5 ) : 
182 -187.[ LUO Ling, WANG Zongming, REN Chunying, et al. 


Models for estimation of grassland production and spatial inversion 
based on MODIS data in Songnen Plain[ J]. Journal of Agricultural 
Engineering ,2010 ,26(5) :182 — 187. | 

苟 其 一 , 董 乙 强 , 安 沙 舟 ,等 .基于 MOD 09GA 数据 的 新 疆 草地 
生长 状况 遥感 监测 研究 [站]. 草 业 学 报 ,2018 ,27(4) :10 -26. 
[XUN Qilei, DONG Yiqiang, AN Shazhou, et al. Monitoring of 


grassland herbage accumulation by remote sensing using MOD 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[32] 


[33] 


[34] 


[35] 


ChinaXiv 合 作 期 刊 


09GA data in Xinjiang[ J]. Acta Prataculturae Sinica, 2018,27 
(4) :10 -26. ] 

青海 省 草原 总 站 . 青海 草地 资源 [M]. 西宁 :青海 人 民 出 版 社 ， 
2012. [ Qinghai General Station of Grassland. Qinghai grassland re- 


sources[ M ]. Xining: Qinghai People’s Publishing House ,2012. | 
金 云 翔 , 徐 斌 , 杨 秀 春 , 等 . 内 蒙古 锡林郭勒 盟 草原 产 草 量 动态 
遥感 估算 [J]. 中 国 科学 :生命 科学 ,2011,41(12):1185 - 
1195. [JIN Yunxiang, XU Bin, YANG Xiuchun, et al. Remote 


sensing dynamic estimation of grass production in Xilinguole , Inner 
Mongolia[ J]. Scientia Sinica ( Vitae) ,2011,41 (12): 1185 - 
1195. | 

ROUSE J W. Monitoring the vernal advancement and retrogradation 
(greenwave effect) of natural vegetation{ R]. NASA/GSFCT Type 
Report, 1974. 

HUETE A,DIDAN K,MIURA T, et al. Overview of the radiometric 
and biophysical performance of the MODIS vegetation indices J]. 
Remote Sensing of Environment ,2002 ,83(1) :195 - 213. 
PEARSON R L, MILLER L D. Remote mapping of standing crop 
biomass for estimation of the productivity of the shortgrass prairie 
[C]//Proceedings of the English International Symposium on Re- 
mote Sensing of Environment, 1972 ,2 ;1375 — 1381. 

ROUJEAN J L,BREON F M. Estimating PAR absorbed by vegeta- 
tion from bidirectional reflectance measurements{ J]. Remote Sens- 
ing of Environment, 1995 ,51(3) :375 —384. 

QI J, CHEHBOUNI A, HUETE A R,et al. A modified soil adjusted 
vegetation index [ J]. Remote Sensing of Environment, 1994, 48 
(2) :119 - 126. 
中 华人 民 共 和 国 农 牧 业 畜 牧 兽 医 司 ,全 国 畜牧 兽医 总 站 . 中 国 
草地 资源 [ M]. 北京 :中 国 科 学 技术 出 版 社 ,1996. [ Depart- 
ment of Agriculture, Animal Husbandry and Animal Husbandry 


and Veterinary of the People’s Republic of China, Ceneral Animal 
Husbandry and Veterinary Station of China. China grassland re- 


sources [ M ]. Beijing: China Science and Technology Press, 


1996. ] 

MKHABELA M S,BULLOCK P,RAJ S, et al. Crop yield forecas- 
ting on the Canadian Prairies using MODIS NDVI data[ J]. Agri- 
cultural and Forest Meteorology ,2011 ,151(3) :385 ~ 393. 
SHRESTHA R,DI L,YU E G,et al. Regression model to estimate 
flood impact on corn yield using MODIS NDVI and USDA cropland 
data layer[ J]. Journal of Integrative Agriculture ,2017 ,16 (2) :98 
-407. 

LOBELL D B, ASNER G P. Cropland distributions from temporal 
unmixing of MODIS data[ J]. Remote Sensing of environment, 
2004 ,93(3) :412 - 422. 

MAKELA H,PEKKARINEN A. Estimation of forest stand volumes 
by Landsat TM imagery and stand-level field-inventory data[ J]. 
Forest Ecology and Management ,2004 ,196(2 -3 ) :245 - 255. 


Ba ESE ET MODIS 的 青海 草地 产 草 量变 化 遥感 分 析 


Grassland yield change in Qinghai Province based on MODIS data 
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Abstract; Qinghai Province is one of the four major pastoral areas in China. It is of great significance to monitor 
the growth of grassland vegetation and accurately estimate grass yields for the sustainable development and ecologi- 
cal protection of the Qinghai pastoral area. Remote sensing estimation of grass yields is mainly based on statistical 
relationship between vegetation index and ground measured data. However, the selection of different vegetation inde- 
xes , various statistical models ,and optional modeling metrics affect the accuracy of the estimation results. Based on 
the ground measured grass yield data of the Qinghai Province and the Moderate Resolution Imaging Spectroradiome- 
ter (MODIS) data,this paper selected six vegetation indexes, five statistical models ,and three modeling metrics to 
test the accuracy of different estimation models. The vegetation indexes included the Normalized Difference Vegeta- 
tion Index (NDVI) , Enhanced Vegetation Index ( EVI) , Ratio Vegetation Index ( RVI) , Difference Vegetation Index 
(DVI) ,Renormalized Difference Vegetation Index (RDVI) , and Modified Soil Adjusted Vegetation Index ( MSA- 
VI). The statistical models consisted of a simple linear model, quadratic polynomial model, power function model, 
exponential function model, and logarithmic function model. The modeling metrics included the annual maximum 
value of vegetation index, cumulative value of vegetation index in growing season, and annual cumulative value of 
vegetation index. First, we explored which vegetation index and estimation model were the most suitable for estima- 
ting the grass yields of different grassland types in the Qinghai Province. Then, we compared the accuracy of esti- 
mated grass yields from the three modeling metrics. Finally ,we selected the optimal grass yield estimation models to 
estimate grass yields of the Qinghai Province in 2015 and 2016 and analyzed grass yield change during 2015 — 
2016. The results showed as follows: (1) among the six vegetation indexes, NDVI was the most suitable vegetation 
index to estimate grass yields in the Qinghai Province. The accuracy of nonlinear models was generally higher than 
that of linear models, especially the exponential models, which are suitable for estimating grass yields of most grass- 
land types. The estimation models based on the annual maximum of NDVI showed the highest coefficient of determi- 
nation ( R°) for almost all grassland types. (2) from the perspective of hay yields, the high-yield area ( >1 200 kg 
- hm~?) was mainly located in the Alpine steppes in eastern Qinghai, and the middle-class area (600 -1 200 kg - 
hm~*) was located in the grass steppes and the Alpine steppes in the north and east of Qinghai. Grasslands with 
low-yield ( <600 kg + hm~?) were located in the Alpine meadows, Alpine steppes, Alpine deserts, and salt mead- 
ows in the western and northern part of Qinghai. (3) the total hay production in the Qinghai Province was 31.60 x 
10° t in 2016, showing 1. 36 % decrease compared to 2015. Among all, the grass steppes and Alpine meadows 
showed the highest yield loss, while in the desert steppes and salt meadows, the yield increased. This study not only 
provided a reference for ecological protection and sustainable development in the pastoral area of Qinghai, but also 
provided research and practices for estimating grassland yield using remote sensing technology. 


Key words; Qinghai Province; grass yield; remote sensing; MODIS; vegetation index 


